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If dark matter (DM) simply consists in a scalar particle interacting dominantly with the Higgs 
boson, the ratio of its annihilation cross section — which is relevant both for the relic abundance and 
indirect detection — and its spin-independent scattering cross section on nuclei depends only on the 
DM mass. It is an intriguing result that, fixing the mass and direct detection rate to fit the annual 
modulation observed by the DAMA experiment, one obtains a relic density in perfect agreement 
with its observed value. In this short letter we update this result and confront the model to the 
recent CoGeNT data, tentatively interpreting the excess of events in the recoil energy spectrum as 
being due to DM. CoGeNT, as DAMA, points toward a light DM candidate, with somewhat different 
(but not necessarily incompatible) masses and cross sections. For the CoGeNT region too, we find 
an intriguing agreement between the scalar DM relic density and direct detection constraints. We 
give the one a region favoured by the CDMS-II events, and our exclusion limits for the XenonlO 
(2009) and XenonlOO data, which, depending on the scintillation efficiency, may exclude CoGeNT 
and DAMA. Assuming CoGeNT and/or DAMA to be due to scalar singlet DM leads to definite 
predictions regarding indirect detection and at colliders. We specifically emphasize the limit on the 
model that might be set by the current Ferrm-LAT data on dwarf galaxies, and the implications for 
the search for the Higgs at the LHC. 

I. INTRODUCTION 

Recently, there has been some effervescence regarding what may be first hints of direct detection of dark matter 
(DM) from the galactic halo. The most recent is related to the CoGeNT experiment, a low background germanium 
crystal based detector, with a rather modest exposure time and detector mass, but very low threshold energy, which 
has announced an anomaly in the form of an excess of events at low recoil energies [T] . Although the collaboration 
clearly leans toward natural radioactivity as the cause for the observed excess, they do put forward the possibility 
that the events may be explained by the elastic collisions of DM from the galactic halo, with a mass in the ^7—10 
GeV range, and a rather large spin independent (SI) cross section on nuclei, a si ~ 7 x 10~ 41 cm 2 . 

Surprizingly, these values for the mass and scattering cross section are not too different from those required to fit 
the DAMA/Libra and DAMA/Nal (DAMA in the sequel) events. DAMA has observed 11 successive cycles of annual 
modulation in the rate of nuclear recoils, with a statistical significance of 8.2<r [2J. These measurements are consistent 
with the signal that would arise from elastic scattering of a WIMP from the galactic halo with the nuclei in the 
detectors, the flux of DM particles being modulated by the periodic motion of the Earth around the Sun 00]. 

There has been much work on the DM interpretation of the recent DAMA data (see e.g. [5-20 ]). In [21] (see also 
[2"2"]h it has been shown that the DAMA results may be explained as being caused by the elastic scattering of a singlet 
scalar DM candidate interacting through the Higgs portal. This particle may be a true singlet scalar [2"3T - f2"5] . or the 
low energy limit of another model 1 . 

II. CONFRONTATION TO COGENT 

In the present short letter, we provide an update on the status of the singlet scalar in the light of the recent data 
from the CoGeNT collaboration, tentatively interpreted as being due to DM. We also include the recent data of the 
CMDS-II collaboration [30 . With 2 events for an expected background of 0.8 ± 0.2 events, the significance of this 
result is low, but it also hints at a rather light DM candidate. 

We adopt the convention of [21 j . and consider just one real singlet scalar S 7 together with a Z2 symmetry, S — > —S, 



A possible implementation is the Inert Doublet Model (IDM), which is another very simple extension of the Standard Model with scalar 

dm [ME!]. 
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to ensure its stability, so that the following renormalizable terms may be added to the SM Lagrangian: 

C3~d*SdpS-~iJ? s !P-^S 4 -\ It H*HS a (1) 
where H is the Standard Model Higgs doublet, and the mass of S is given by 

m%=ii% + Alv 2 (2) 

with v = 246 GeV. Both annihilation into SM particles and scattering with nuclei are through the coupling to the 
Higgs particle h, respectively in the s and t-channel. Annihilation through the Higgs is S-wave, and the cross section 
for scattering on nuclei is purely of the SI type. For a DM candidate lighter than the Higgs, 7715 <C mh, the ratio of 
the annihilation and scattering cross section depends only on m<j, 

v a{SS -> gKej = n - m \ ( m s ~ m W 2 ( ) 

Y <j(SN -> SN) Y f 2m2 N^r m s 

where n c = 3(1) for quarks (leptons), and \i r = msmN/(iTis + t^n) is the nucleon-DM reduced mass. The factor / 
parametrizes the Higgs to nucleons coupling, fniN = (N\ J2q TO g9<7l A) — gtiNN^, and we consider 0.2 < f < 0.4 (see 

Equation ^ shows that the mass of the DM candidate is fixed for a given relic abundance and SI scattering cross 
section [211 I23H25] . In turn, direct detection experiments may determine both the SI cross section and the mass of the 
DM, modulo the astrophysical uncertainties regarding the local density and velocity distribution of the DM. A priori 
there is little chance that these constraints may be met by singlet scalar DM, but as Figure [T reveals, the model may 
be in agreement with CoGeNT — which is the main result of this letter — or, as shown in [5TJ[55], with DAMA. Since 
there is a gap between the CoGeNT and DAMA (with channelling) regions, that scalar DM agrees with CoGeNT 
does not trivially derive from the fact that it may agree with DAMA. We emphasize that this result, as for DAMA, is 
specific to a scalar particle with scalar couplings to SM fermions. For instance, annihilation through the Higgs portal 
would be P-wave suppressed for a fcrmionic singlet DM candidate, and other interactions, as is the case for a light 
neutralino (32], are necessary to agree with the direct detection data (see also e.g. [33"rf36] ) . 

The gap between CoGeNT and DAMA may be reduced, either assuming that channelling is less effective than 
what is advocated by the DAMA collaboration (which has the effect of raising the DAMA region — but not reducing 
the tension with exclusion limits), or assuming that the CoGeNT excess is partially contaminated by some natural 
radioactivity (lowering the CoGeNT region) or a mixture of both. One may also adjust the properties of the halo or 
the DAMA spectral data [33], but we have refrained from doing so. 

Also, both regions are excluded by the most stringent limit set by XenonlO (the dotted line in Figure [l]). However, if 
the scintillation efficiency — the measure of the fraction of energy from the recoiling nuclei that goes into scintillation 
light — is actually lower than that used by the collaboration (the long-dashed line in Figure [l]), as advocated in [40] . 
there is a region of CoGeNT which may be consistent with all experimental constraints (and the singlet scalar DM 
candidate). 

In the scenario we consider, the coupling A^ to the Higgs must be fairly large to explain both the relic abundance 2 
and the direct detection data. For mh = 120 GeV (m^ = 180 GeV), one typically requires A^ ~ —0.2 (resp. ~ —0.45) 
for a DM candidate of mass mg ~ 8 GeV. For the same choice of parameters, fis ~ 100 GeV, which implies fine tuning 
of the parameters at the level of the percent, which is not unbearable in our opinion. Also there is no mechanism to 
naturally stabilize the mass of the scalar at a scale of few GeV. Note however that, if neither CoGeNT nor DAMA 
could be fitted with /i| ~ at tree level, in the Inert Doublet Model at one loop, the DM mass and coupling ranges 
required by CoGeNT and/or DAMA may be compatible with dynamical electroweak symmetry breaking induced by 
the inert doublet [4"3"] . 

A light scalar dark matter candidate coupled to the Higgs has potentially many other signatures or implications, a 
large flux of gamma rays from dark matter annihilations [101 121] , a large flux of neutrinos from capture by the Sun, 
which may be constrained by Super-Kamiokande FI4TI41)] , or anti-protons and anti-deuterons in cosmic rays [5] H7] . 



2 In the calculation of the relic abundance, we neglect the possible effect of the QCD phase transition, which is supposed to have taken 
place around ~ 150 MeV (see, for instance, |41|). Since typically Xf = mg/Tf ~ 20, for mg ~ 6 GeV, for instance, the freeze-out 
temperature is Tj ~ 300 MeV 2> 150 MeV and thus the QCD phase transition is irrelevant for the range of mass we consider. Otherwise, 
the QCD phase transition might increase the relic abundance by at most a factor of 0(2) (see for instance |42|). Such an increase could 
be compensated by an increase of by a factor of ~ 1.4 and this, in turn, would require a decrease of the parameter / by the same 
factor. 
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Figure 1: SI cross section (a„) vs scalar singlet mass (ms), for pom = 0.3 GeV/cm 3 and a standard Maxwellian velocity 
distribution (with mean velocity 220 km/s and escape velocity v esc = 650 km/s, see our conventions in The green region 

corresponds to CoGeNT (minimum x 2 , with contours at 90 and 99.9% C.L.), for which we have assumed that the excess at low 
recoil energies is entirely due to DM (assuming a constant background contamination). The DAM A regions (goodness-of-fit, 
also at 90 and 99.9% C.L.) are given both with (purple/orange) and without (purple, no fill) channelling. The blue region 
corresponds to the CDMS-II two events, at la, which we obtained following the procedure of |37j . The blue (short-dashed) line 
is the 90% C.L. exclusion limit from CDMS-Si [35]. The black dotted line is the 90% C.L. exclusion limit from the XenonlO 
2009 data set, using their scintillation efficiency [39|, as also considered in [37]. The long-dashed line is based on the same 
data but using instead the smaller scintillation efficiency advocated in [3D] (central value, at la the corresponding exclusion 
can be found in 34 ). Finally, the brown lines (continuous) encompass the region predicted by the singlet scalar DM model 
corresponding to the WMAP range 0.094 < Q, DM h 2 < 0.129, for 0.2 </ < 0.4. 
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10 GeV 
BR(SS -> t+t-) ~ 10% 
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6 GeV 
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BR(SS -^bb + cc) ~ 80% 
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1.0 x 10~ 9 


2.8 x 10~ 9 


1.7 x 10~ 9 



Draco 



Table I: Comparison between the expected gamma-ray flux from a light scalar and the 95%C.L. limits given by the Ferrm-Lat 
collaboration, Figure 2 in [48]. For the 10 GeV candidate the limits are extracted assuming annihilation into bb with a BR of 
100%. The limits for the 6 GeV candidate are our extrapolations, assuming BR=80% BR in bb and BR=20% in t + t~. 



Following a suggestion made in |34j . we may confront the model to data on the gamma flux from dwarf galaxies 
recently released by the Fermi-L AT collaboration @8J. The analysis in [?S] gives, for various dwarf galaxies, the 95 
% C.L. limit on the total flux $ of gamma rays (with energy between 100 MeV and 50 GeV) that may be produced 
through annihilation of dark matter. The published analysis, which is quite sophisticated, is limited to candidates 
with a mass larger than 10 GeV. However the spectrum of photons is quite similar for slightly lighter candidates 
(see Figure |2|, so we expect the constraints to extrapolate smoothly for, say, a 6 GeV candidate. For the sake of 
illustration, we consider the limits from two representative dwarf galaxies, Draco and Ursa Minor [IS]. In Table I, we 
give the predictions for the singlet scalar model for candidates with mass 6 and 10 GeV, assuming the NFW profiles 
as used by the collaboration (see Table 4 of @8]), and for av = 2.5 • 10~ 26 cm 3 -s _1 . 
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Figure 2: EdN/dE spectra produced by Pythia8.1 for a scalar singlet of mass 6 GeV (blue, long dashed), 8 GeV (dark blue, 
short dashed) and 10 GeV (black, continuous line). 

The 10 GeV (6 GeV) scalar interacting through the Higgs has a BR into t + t~ of order 10% (20%). For the limits 
on the flux, we refer to Figure 2 in [48], which display the limits obtained for a BR=100% into bb, and for BR=80% 
into bb and BR=20% into t + t~ . For the 10 GeV candidate, it is a good approximation to take the limits assuming a 
BR of 100% in bb (Table I). As shown in Figure 2 in g5], the BR~ 10% in t+t~ should give a limit on the flux that 
is O(10%) stronger, at most. For a 6 GeV candidate, we use the figure with BR~ 20% in t + t~ and BR~ 80% in bb 
and naively extrapolate the curve down to 6 GeV. In all cases the predicted fluxes are larger than the limits. From 
the numbers given in Table I, we thus tentatively conclude that the scalar singlet model with a NFW profile may be 
excluded at 95 % C.L. by data from dwarf galaxies, and that it would be interesting to extend the analysis of [H] to 
lighter WIMP candidates. 

To close the discussion on indirect constraints, we mention the fact that the annihilation and mass of the scalar 
singlet candidate makes it also a very natural candidate to solve the primordial 6 Li problem (see Figure 3 in [4T?]). 

Last but not least, a light WIMP in the form of a scalar coupled to the Higgs would imply that the Higgs mostly 
decays into a pair of dark matter particles [21] [24] [25] [50] . We would like to point out that the DAMA and CoGeNT 
regions could be distinguished from a measurement of the invisible Higgs decay branching ratio. For to/, = 180 GeV 
the effect is striking: as Figure [3] shows, taking into account the CDMS-Si limit and the WMAP region, the DAMA 
region gives 60% <BR< 70%, while for the CoGeNT region, one has 75% < BR < 90%. This difference is larger 
than the expected ~ 10% LHC sensitivity on the invisible branching ratio [BT] . For rrih = 120 GeV, the difference 
is much reduced, because the invisible channel largely dominates the decay width: we get 98% < BR < 99% and 
BR > 99% for DAMA and CoGeNT respectively. 

III. XENONIOO EXCLUSION LIMITS 

Shortly after CoGeNT, the XenonlOO collaboration has released an analysis of their first data, collected over only 
11.2 days [52 . The collaboration has found no event consistent with dark matter, and it has the world best SI 
exclusion limits for dark matter mass ^ 80 GeV (Figure 5 in |52|). More surprizing, at first sight, is the limit set 
on lighter WIMPs, with mass % 10 GeV. In particular the CoGeNT and DAMA (with and without channelling) are 
excluded with 90% confidence level 3 . In this section, we present our own analysis of the low mass region, following 



3 That something non-standard is being done in this region may be appreciated by comparing Figure 5 of [52] to the preliminary results 
presented at recent conferences (see for instance 53 ) and is subject to discussion |54W56| . 
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the guidelines provided in Refs. |521 154H56] . to which we refer for more details. For light dark matter candidates, the 
relevant features of the analysis of XenonlOO are as follows. 

To discriminate dark matter (and neutron) collisions with Xe from electron and gamma events (a large source of 
background) the experiment relies on photoelectrons (PE) produced by scintillation (SI signal). The relation between 
the mean number of PE, npE(E nr ), and the recoil energy E nr (in keV) is given by the so-called scintillation efficiency 
(Leff). Which Leff to use, and how to extrapolate Leff at low E nr (where no measurements of Leff exist) is an 
important issue [54] . The experimental situation on Leff is summarized in Figure 1 of Ref. [55] . The exclusion limit set 
by XenonlOO in [53] is based on the best fit to current experimental data (LeffMed here), which gives Leff s» 0.12 at 
small nuclei recoil energies E nr . Furthermore Leff is set to zero for E nr % 1 keV. A more conservative choice (LeffMin 
here, corresponding to a lower 90% C.L. fit to the data) gives a Leff which decreases monotonically with E nr and 
vanishes at E nr < 1 keV. The Zeplin experiment (also a Xe experiment) uses a different Leff, which is essentially zero 
below 6-7 keV (LeffZep here) [57]. 

In Figure |4j following [55], we show the theoretical event rate (per kg and per day) for a mass = 10 GeV and 
(Tsi = 10~ 41 cm 2 candidate, typical of the DAMA region. We use v esc = 544 km-s -1 as the XenonlOO collaboration. 
In the same Figure, the size of the bins correspond to OPE, 1PE, etc., given here for LeffMin. The threshold depends 
on the acceptance of the detector. For 3PE the acceptance is about 50%, and about 70% for 4PE (Figure 3 of 
Ref. [55]). A standard analysis would give no event with 3PE or more. However, because of fluctuations (assumed 
to be poissonian) in the number of PE produced, and because the event rate is exponentially rising at low E nr (for 
elastic scattering), the rate is much larger at 3PE and 4PE, as shown by the histograms in Figure |1J obtained by 
convolution of the theoretical rate with a Poisson distribution of mean npE(E nr ). For the sake of comparison we give 
the result of our calculation (black histogram), and that of XenonlOO (green histogram), see Figure in [55] . We have 
slightly more events at 3PE (hence we will have slightly stronger limits), but otherwise the agreement is good. 

Our exclusion limits are shown in Figure [5] where, in the left panel, we give the limits for a threshold at 3PE 
(consistent with Xenonl00|55J) and in the right panel the one for the more conservative choice of 4PE (consistent 
with Ref. [54 ). From left to right, the exclusion limits correspond to LeffMed (short dashed), LeffMin (continuous) 
and LeffZep (long dashed). For the latter, the cutoff in Leff is such that the effect of Poisson fluctuations from events 
at low E nr is negligible, and one recovers the limit from a standard analysis. 

An important feature we observe is that a large fraction of the events at 3PE/4PE are associated to Poisson 
fluctuations from the small E nr region, where the rate is very large, but less than 1PE is expected on average. As 
a result, if one increases the E nr at which Leff vanishes, the number of low energy events that are lost (and which 
consequently do not contribute to higher energy bins) quickly increases 4 . This effect is illustrated in the right panel of 
Figure |4j where we show the number of events above 3PE (resp. 4PE) for candidates with a si = 10~ 40 cm 2 , both for 
LeffMed and LeffMin, but changing the cutoff E nr at which they are assumed to vanish. The effect is most sensitive 



4 This loss is the reason why the number of events in the first bin is about a factor 2 smaller than the theoretical number of events. 
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Figure 4: Left panel: Effect of Poisson fluctuations on the expected signal (continuous curve) from a candidate with mass = 10 
GeV and a si = 10~ 41 cm 2 , assuming LeffMin (in green the XenonlOO bins, in black our bins - see text). The two vertical 
dashed lines correspond to the 3PE and 4PE thresholds. 

Right panel: Effects of changing the cutoff in the recoil energy (see text) on the number of events above 3PE (in black) and 
4PE (in red) for LeffMed (short dashed), LeffMin (continuous), for a candidate of 6.5 GeV for LeffMed and of 7.5 GeV for 
LeffMin, both for a si = 10 -40 cm 2 . The horizontal dashed line corresponds to the 90% C.L. exclusion limit, which corresponds 
to 2.3 events according to Poisson statistics. 




m s (GeV) m s (GeV) 

Figure 5: XenonlOO exclusion limits with 90% C.L., with threshold at 3 PE (in black, left panel) or 4 PE (in red, right panel). 
The curves correspond respectively to the LeffMed (short dashed), LeffMin (continuous) and LeffZep (long dashed) scintillation 
efficiency — see text. For the sake of comparison, we have taken v esc = 544 km-s -1 like the XenonlOO collaboration. The 
blue (dot-dashed) lines correspond to our predicted exclusion limit for XenonlOO, using LeffMin and for an exposure of about 
1 ton-days, assuming zero event. 

for LeffMed, where for a cutoff at E nr = I keV the candidate is excluded at 90% C.L., while it is allowed if the cutoff 
is at E nr — I .5 keV. In our exclusion limits we use the function LeffMed and LeffMin with a cutoff at E nr = 1 keV. 
In our opinion, the method followed by the XenonlOO collaboration is sound, but is very sensitive (as discussed in 
to the choice of Leff at low E nr . Here we would like to emphasize that not only the shape of Leff, but also the 
E nr at which Leff is cutoff, is a critical issue in setting exclusion limits. This being said, one must admit that there 
is a tension between the CoGeNT (and a fortiori DAMA) regions and the limits set by the XenonlOO experiment, 
and this with only limited exposure. Assuming that no events are seen in the data, we give in Figure [5] our predicted 
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exclusion limit of XenonlOO, using LeffMin, for an exposure corresponding to 1 ton-days, as reported in [S3] . 
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